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20.  Abstract  -  Contd.  ~  —  ■— 

This  report  describes  the  development  of  a  model  that  may  be  used  to 
estimate  the  parameters  of  precipitable  rain  distributions  from  inputs  of  liquid- 
water-content  and/or  measurements  of  radar  reflectivity  coupled  with  stand¬ 
ardized  cloud  physics  relationships. 

The  model  is  developed  from  an  equation  set  based  on  an  exponential  dis¬ 
tribution  function.  Aircraft -acquired  data  are  used  to  verify  the  conformity  of 
rain  distributions  to  exponential  shapes.  Empirical  relationships  provided  by 
these  data  verified  the  existence  of  a  nondimensional  predictable  entity  (A Dm  = 
distribution  slope  times  maximum  drop  size),  which  provides  improved  esti¬ 
mates  of  rain  distributions  from  predicted  or  measured  values  of  liquid-water- 
content  and  radar  reflectivity. 

Tables  listing  the  variations  in  the  size  distributions  during  three  rain 
situations  are  given  in  Appendix  A. 
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A  Model  for  the  Estimation  of  Rain  Distributions 


I  IM'KOIH  C’.TION 

Although  detailed  i-j-iit  rophysical  hvdrometeor  definitions  an>  not  a  necessary 
nyutn'iuf'it  in  tin-  initial  development  of  dynamic-  models  for  future  application 
it.  meteorological  f-  e.cast  ing.  a  n  \  effort  designed  to  predict  or  describe  situa¬ 
tion.-  in  which  pree  ipir -,t  ion  is  present  has  to  be  able  to  produce  realistic  facsim¬ 
iles  of  hydrometeor  type,  concent  rat  ion,  and  size  distribution.  This  is  particu¬ 
larly  nei-.  -s’iiy  in  fore,  testing  scenar  ios  where  a  choice  has  to  be  made  between 
various  systems  'hat  have  weather-dependent  operational  efficiencies,  such  as 
in  aircraft  safety,  optica!  guidance  (smart  weapons),  and  communications.  Thus, 
it  is  extremely  important,  at  'his  time,  before  our  dynamic  models  reach  the 
point  of  be.  onting  operational,  to  develop  the  methods  by  which  information  of 
tins  Hue  cm  be  included  in  future  large-scale  models  so  that  they  may  be  used 
as  effective  forecasting  tools. 

This  report  describes  the  development  of  a  mathematical  model  or  equation 
set  that,  ran  reasonable  describe  the  parameter  s  of  pi  ecipltable  liquid  hydrome- 
teors  for  the  subsequent  inclusion  in  more  comprehensive,  large-scale  models. 
The  effort  i  neompassed  the  analysis  of  aircraft  distribution  data  for  the  empirical 
determination  of  mathematical  relationships  to  describe  the  sizes  and  shapes  of 

(He.-eived  for  publication  27  January  1988) 
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hvdrometeor  populations  and  the  documentation  of  the  variabilities  of  rain  distri¬ 
bution  parameters  that  may  be  used  to  establish  probabilities  of  occurrences. 


2.  DERIV  ATION  OF  EQUATIONS 

1  9 

It  has  been  demonstrated  by  Marshall  and  Palmer,  Marshall  and  Gunn,  “ 
Imai  et  al,  Gunn  and  Marshall,  ^  Ohtake  and  Henmi0  and  others  that  the  size 
distribution  properties  of  raindrops,  snowflakes,  and  ice  crystals  of  preripitable 
size  can  be  reasonably  described  by  a  distribution  function  of  exponential  type. 
This  distribution  function  specifies  that  the  number  concentration  of  the  hydrome¬ 
teor  par  ticles  will  decrease  with  increasing  diameter-  (or  equivalent -melted 
diameter)  in  the  manner 


N  -  N“ e'AD  ,  (d  <  D  <  D  ) 
0  m 


No. 


-:i 

m 


-1 

mm 


(1) 


where  Nq  is  the  number  per  cubic  meter  per  millimeter  bandwidth  at  the  zero 
intercept  of  the  semi-logarithmic,  number-density  plot,  A  is  number  per  milli¬ 
meter  bandwidth  and  is  the  slope  of  the  number-density  dist  ribution,  and  D  is  the 
drop  diameter  in  millimeters. 

The  equation,  as  applied  herein,  is  presumed  to  be  descriptive  only  between 
the  truncation  limits  D  d  (a  minimum  diameter)  and  D  =  D  (a  maximum  diam¬ 
eter).  This  subject  of  the  double  truncation  of  an  exponential  distribution  function 
has  been  previously  discussed  by  Sekhon  and  Srirastava.  l> 

The  total  number  of  hydrometeors  (N,p)  in  a  population  described  by  Kq.  (1)  is 


N.„ 


N  dl) 


No.  m 


(2) 


1.  Marshall,  J.  S.  ,  and  Palmer,  W.McK.  ( 1 P 4 8)  The  distribution  of  raindrops 

with  size,  J.  Meteorol.  5:  ltia- 1  (»(». 

2.  Marshall,  J.  S.  ,  and  Gunn.  K.  1,,S.  (1952)  Measurement  of  snow  parameters 

by  radar,  J.  Meteoiol.  9:222. 

3.  Imai,  I.,  Fujiwara,  M.  ,  Ichimura,  I.,  and  Toyama,  V.  (1955)  Radar  reflec¬ 

tivity  of  falling  snow,  Pap.  in  Meteorol.  and  Geophys.  (Japan)  (1:120-139. 

4.  Gunn,  K.  I,.  S. ,  and  Marshall,  J.S.  (1958)  The  distribution  with  size  of  aggre¬ 

gate  snowflakes,  J.  Meteorol.  15:452(479). 

5.  Ohtake,  T. ,  and  Hennri,  T.  (1970)  Radar  reflectivity  of  aggregated  snow  - 

“■lakes.  Preprints  of  papers  presented  at  the  14th  Radar-  Meteorology 
C'onfm-ence,  Tucson,  Arizona,  17-20  November  1970,  pp.  209-211. 

(i.  Sekhon,  R.S.  ,  and  Srivastava,  R.C.  (1970)  Snow  size  spectra  and  radar 
reflectivity,  J.  Atmos.  Sc i .  27:299-307. 
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T  A 
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where  is  a  "truncation  ratio"  specified  by 


r 


N  dD 


(4) 


/ 


N  dD 


which  becomes 


1  \  -  ° 


-<JA 


■DmA 


(5) 


The  liquid-water-content  or  mass  (M)  of  the  hydrometeor  populations  de¬ 
scribed  by  Kq.  (1)  is  distributed  with  diameter  as  a  function  of  the  third  moment 
of  Kq.  (1),  or  as 


. ,  it  , . -3  »3  -AD  -3  -1 

Rip,  -  nr  <  10  p  N„  O  e  p  m  mm 
On  r  ul  0  ” 


(«) 


where  d<  0<  Dm  and  p  is  the  density  of  liquid  water-  in  g  cm 
The  total  liquid-water-content  of  the  population  is 

D 


m7 


M  p  d  D  g  m 


(7) 


which,  from  b'q.  (6)  and  integration,  yields 


tr  •  10*3  p  N  r(4)r  „ 

u  0  M  -3 

1V1  = - j -  g  nn 

BA  4 


(8) 


where  F  (4)  is  the  gamma  function  of  4  and  r^  is  a  truncation  ratio  for  liquid- 
watei  -c  ontent  given  by 
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rM  =  \  )e‘dA  [(dA)3  +  3(dA)2  +  6  dA  +  6]  -  e  °mA  * '3 


l(DmA)J 


+  3(DmA)  +  6DmA  +  6]  | 


(10) 


The  distributed  values  of  the  radar  reflectivity  factor  (Z)  for  the  hydrometeor 
populations  described  by  Eq.  (1)  are  specified  by 


Z„  =  N-  D6  e‘AD  .  (d  <  D  <  D  )  ,  mm6  m'3  mm-1  . 

DO  •  m 


(11) 


The  total  value  of  the  radar  reflectivity  factor,  for  the  entire  population  of 
hydrometeors,  is 


D 


f  *0 


dD  mm^  m  3 


(12) 


or,  from  Eq.  (11)  on  integration, 


Z  = 


N0r(7)  rZ  6 


(13) 


A 


where  T  (7)  is  the  gamma  function  of  7  and  r^  is  the  truncation  ratio  for  the  radar 
reflectivity  factor  as  defined  by 
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It.  is  seen  th.it  I)  is  nons'-pn'  iibh-  in  this  equation.  However,  the  equation 

he  ‘Xt  li!v  solve  1  it'-  ’  :  i;,l  -  in  >  or  or.ee  information  about  d,  U  and  A 

m 

v-ilahl.  . 
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The  parameter  N^  may  be  eliminated  between  Eqs.  (8)  and  (13)  to  provide  an 
expression  for  the  "exponential  slope"  of  the  distribution  function  of  Eq.  (1). 
Thus,  after  evaluation  of  p  ,  T  (4),  and  T(7) 

Iw.' 


/M 

1  3 

-1 

U 

mm 

\Vr  can  also  eliminate  the  parameter  A  between  Eqs.  (8)  and  (13)  to  obtain 
an  expression  for  should  we  wish. 

Equations  (1)  through  (20)  constitute  a  descriptive  equation  set,  or  model, 
that,  with  several  assumptions  as  will  be  discussed  here,  can  be  solved  in  closed 

form. 

I  he  distribution  equations,  Eqs.  (1),  (r,),  and  (11),  for  a  truncated  model  are 
identical  to  those  for  a  nontruncated  model  except  for  the  recognition  that  the 
truncated  equations  have  significance  onlv  between  the  diameter  limits  D  -  d  to 
I)  1)  . 

The  totals  equations,  Eqs.  (3),  (8),  and  (13),  differ  from  those  of  a  nontrun¬ 
cated  model  in  that  the  former  contain  the  truncation  ratios,  r.,,  r„,,  and  r_  . 
ns  defined  by  Eqs.  (5),  (10),  and  ( 1  ">).  These  truncation  ratios  are  seen  to  be 
functions  of  1,  I)  ,  and  A. 

The  equations  presented  in  this  section,  as  previously  mentioned,  are  perti¬ 
nent  to  precipitable  liquid  hydrometeors  onlv,  namely  rain  or  the  result  ing  melted 
drops  from  snow /ice  particles.  This  particular  investigation,  henceforth,  will 
only  be  concerned  with  tain.  It  is  planned  to  continue  these  studies  into  the 
snow/ ice  region. 


3.  DATA  ANALYSIS 

P’e  EY  (  archives  were  searched  to  find  the  best  rain  situation  to  analyze  in 
or  le:  to  investigate  the  validity  of  the  preceding  equation  Set.  The  case  selected 
was  on'  of  widespread  rain  that  was  sampled  by  a  MC-130K  inst -umented  aircraft 
nett  ialla  iega,  Alabama  on  28  I  "brur  tv  11*77,  where  measurements  of  precipita¬ 
tion  we*,  taken  continuously  for  24  min,  from  2222  to  2248  GMT.  Data  used  in 
this  analvsis  were  from  a  FMS  1-U*  precipitation  probe  that  counted  raindrops 
I  tom  0.2-  to  4.  til*  -  mm  ilia  met  e  r  and  classified  them  by  size  in  13  channels  of 


Knol lenberg,  11. G.  (lf*70)  Fhe  optical  array:  an  alternative  to  scattering  or 
extinction  for  airborne  particle  size  determination,  J.  Appl.  Meteor.  9  (No. 
11:80-103.  - - — — - 
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—  0.3-mni  width.  I'he  mpasurcmpnts  ivpro  takon  at  a  1-sec  time  resolution. 

This  particular  rase  had  data  recorded  in  each  of  the  1440  1 -second  samples. 

This  rain  situation  has  been  previously  analyzed  using  a  nondimensional 

3  > 

technique  (Plank,  Berthel,  and  Delgado)  to  test  the  <  nontial  distribution  as¬ 
sumption  for  tain.  That  investigation  determined  the*  distributions  of  liquid  hy- 
dtometeors  could  be  described  by  exponential  functions  with  the  distributions 
becoming  more  exponential  in  form  as  averaging  time  increased. 

Similar  findings  are  illustrated  in  Figure  1.  The  number- density  distribution 
plots  for  the  initial  1,  a,  10,  a0,  100,  rind  aOO  sec  of  the  sampling  period  show 
the  number  of  hops  (normalized  to  a  cubic  meter  pet  millimeter  bandwidth)  in 
each  PMS  1-D  hanm-l  plotted  vs  the  midpoint  linmetet  of  each  class.  The  solid 
lines  are’  best  fit’  litres  h-rived  from  least  -squat  >■>  regression  analysis.  The 
dashed  litres  am-  slntie;  eaU  ulate  l  using  Fq.  (20).  (Determination  or  the  slope  in 
•Iris  mature:  is  Its  uss.-l  in  more  detail  later  in  this  section.  *  The  points  in  th< 
firs:  plot  are  from  basic  t-se.  lata,  f!"  rthers  are-  mean  values  found  Ir  sum¬ 
ming  the  t . •suits  of  each  channel  for  »he  basic  period  and  lividing  In  the  averaging 
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Figure  2.  Plot  of  the  Mean,  Minimum,  and  Maximum 
Values  of  \.j,  vs  Averaging  Time  on  2d  February  1977 
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Figure  3.  Plot  of  the  Mean,  Minimum,  and  Maximum 
Values  of  D  vs  Averaging  Time  on  23  February  1977 
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and  the  number  of  drops  per  cubic  meter  allows  this  value  to  be  calculated,  as¬ 
suming  each  drops  to  be  spherical,  as 


M;  ~  <  10'3  pw  13?  X.  g  m*3  ,  (21) 

where  "  i"  is  the  class  designation,  p  ,  is,  once  again,  the  density  of  liquid  water 
-3  w 

(g  cm  ),  D.  is  the  middiameter  of  class  "i",  and  Nj  is  the  number  of  drops  per 
cubic  meter  in  the  class. 

The  total  M  of  a  particular  distribution  is  found  by  summing  the  contents  of 

each  class  as 


i  - 15 

M  -  j-  <  10*3  pw  ^  D3  N  g  m'3  .  (22) 

i  =  1 

The  equivalent  radar  values  can  be  calculated  in  similar  fashion  as 

Z.  =  d!5  N.  mm6  m*3  (23) 


and 


i=  15 

Z  =  D?  N.  mm6  m~3  .  (24) 

i=l 

The  means,  minimums,  maximums,  and  standard  deviations  for  M  and  Z  are 
plotted  in  Figures  4  and  5  for  various  averaging  periods. 

Of  these  four  parameters,  D  displays  the  least  amount  of  variability  (Table 
A3)  within  this  specific  case  example.  However,  all  will  differ  widely  in  different 
rain  situations.  Thus,  none  can  be  considered  a  predictable  quantity. 

Examination  of  the  equations  in  Section  2  reveal  other  parameters;  Nq  ,  Dq  , 
and  A  and  the  truncation  ratios  r^ ,  r^,  and  r^,  that  can  be  termed  "derived 
variables"  in  that  they  are  derived  from  calculations  using  the  data  from  the  PMS 
1-D  and  the  assumption  of  exponential  shape.  The  single,  most -identifying  feature 
of  any  exponential  distribution  is  A,  as  it  is  an  integral  part  of  all  the  equations 
cited  in  Section  2.  This,  the  slope  of  the  number -density  distribution,  can  be 
determined  through  knowledge  of  M  and  Z  and  utilization  of  Eq.  (20)  as  previously 
mentioned.  The  slope  can  also  be  derived  by  using  N^,  with  M  or  Z  by  eliminating 
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N0  and  equating  Eqs.  (1)  and  (8)  or  Eqs.  (1)  and  (13).  Since  this  investigation  is 
primarily  concerned  with  the  prediction  of  realistic  rain  distributions  from  fore¬ 
casts  of  dynamic  models,  we  chose  to  use  the  M  and  Z  as  these  parameters  are 
the  more-likely  candidates  for  prediction  modeling  (this  subject  will  be  discussed 
further  in  Section  4). 

The  dashed  lines  in  Figure  1  show  the  A's  that  were  calculated  using  Eq.  (20) 
for  each  of  the  distributions.  This  method  of  determining  A  assumes  a  full  dis¬ 
tribution  or.  in  other  words,  a  representative  number  of  drops  in  each  class  that 
will  fulfill  the  requirements  imposed  by  the  M  and  Z  parameters.  The  least- 
squares  method,  which  is  dependent  on  distributed  numbers,  does  not  take  into 
account  the  dependency  of  M  on  diameter  to  the  third  power  [Eq.  (21))  and  Z  on 

diameter  to  the  sixth  power  [  Eq.  (23)].  The  small  averaging  intervals  represent 

3 

extremely  small  sampling  volumes  (<~  1  m  )  and  some  of  the  1-D  classes  are 

9  ... 

deficient  or  devoid  of  drops  (Plank  and  Berthel).  Therefore,  it  is  not  surprising 
that  the  two  methods,  Eq.  (20)  and  the  least  squares,  differ  at  the  small  sampling 
intervals  but  tend  to  agree  more  as  the  averaging  period  becomes  larger.  (Simi¬ 
lar  results  are  obtained  when  using  Nj,  in  the  calculation  of  A. ) 

Figure  6  shows  the  means,  minimums,  maximums,  and  standard  deviations 
that  were  calculated  using  M  and  Z  in  Eq.  (20)  for  different  averaging  periods. 

The  calculations  of  A  in  Figure  6  used  the  experimentally  determined  values 
of  Dm  in  solving  for  rM  and  r^  [  Eqs.  (10)  and  (1 5)1 .  In  any  scenario  using  pre¬ 
dicted  values  of  M  and  Z,  Dm  will  be  an  unknown  quantity.  Thus,  some  judgement 
or  assumption  has  to  be  made  as  to  the  largest  drop  size  that  would  most  likely 
be  present.  When  the  equations  for  the  truncation  ratios  of  M  and  Z  are  consid¬ 
ered  [Eqs.  (10)  and  (15)1,  it  is  apparent  that  the  two  unknown  values  are  A  and  D[r) 
since  d  can  be  defined  as  the  diameter  where  a  drop  is  of  sufficient  size  to  become 
precipitable.  (A  precise  definition  of  d  is  not  always  necessary  since  small 
changes  do  not  significantly  effect  M,  Z,  or  the  calculated  A •  But,  it  does  have  a 
decided  effect  in  the  calculation  of  the  total  number  of  drops  because  of  the  nega¬ 
tive  slope  of  the  exponential  number  distribution.) 

Multiplication  of  the  A  and  Dm  parameters  forms  a  nondimens ional  entity 
(AD  )  and,  if  d  and  A  Dm  are  assigned  values,  leaves  a  single  unknown,  A,  in 
Eqs.  (10)  and  (15).  This  manipulation  now  allows  A  to  be  determined  through  the 
solving  of  Eq.  (20)  by  the  trial-and-error  method.  Figure  7  shows  the  ADm 
means,  minimums,  maximums,  and  standard  deviations  for  this  sample  case  over 
different  averaging  periods. 


9.  Plank,  V.G.,  and  Berthel,  R.O.  (1982)  A  descriptive  double-truncated  expo¬ 
nential  model  for  hydrometeors  of  precipitable  size,  Preprints  of  papers 
presented  at  the  Conference  on  Cloud  Physics.  Chicago,  Illinois,  15-18 
November  1982,  pp.  190-194,  AFGE-TR-82-0347,  ADA122038. 
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Figure  6.  Plot  of  the  Mean,  Minimum,  and  Maximum 
Values  of  A  vs  Averaging  Time  on  23  February  1977 


Figure  7.  Plot  of  the  Mean,  Minimum,  and  Maximum  Values  of  AD 
vs  Averaging  Time  and  Sampling  Volump  on  23  February  1977 


All  the  parameters  discussed  in  this  section  vary  considerably  within  any 
specific  rain  situation  with  the  variance  being  dependent  upon  averaging  period. 
These  same  parameters,  with  the  exception  of  ADm,  can  show  substantial  varia¬ 
tions  between  different  rain  cases.  AD|T],  although  displaying  a  dependence  upon 
averaging  period,  tends  to  have  the  same  approximate  value  in  all  rain  situations. 
Thus,  it  can  be  used  as  a  predictable  quantity. 


4.  DISCUSSION 

Our  presumptions  about  the  parameters  A  and  Dm  were  the  following.  For 
any  single  hydrometeor  sample,  irrespective  of  spectral  type,  that  is,  exponen¬ 
tial,  bi-modal,  etc.,  it  has  been  demonstrated  by  Plank,  Berthel,  and  Barnes^ 
that 


M  =  k  Z'  5 


g  m 


(25) 


where  k  is  the  so  called  "kappa  factor"  that  can  be  evaluated  from  aircraft  data. 
When  this  equation  is  substituted  into  Eq.  (20). 


A 


61. 2 


1/3 


mm 


(26) 


For  a  family  of  hydrometeor  samples,  we  presumed  that  an  M  vs  Z  relation 
existed  that  was  of  power  function  form, 


IV1  =  a  F.13  g  m 


(27) 


and  that  was  assumed  to  be  known  either  from  literature  information  or  from 
regression  analyses  performed  on  .ioint  aire -aft -radar  data.  If  Eq.  (27)  is  sub¬ 
stituted  into  Eq.  (20), 


A 


61.2 


1 3 b  ,.h-l.'3h 
a  i\l 


mm 


(76) 


10.  Plank,  V.G.,  Berthel,  K.O.,  and  Barnes,  Jr.,  A.  A.  (1080)  An  improved 
method  for  obtaining  water  content  values  of  ice  hvdi  omet enrs  from  air¬ 
man  and  radar  data,  J.  Appl.  Met  enrol.  1 0: 1 203  - 1200,  A  EG  I  ,-Tlt-il  1 -001 1 
AD  AOP4328.  '  ~  .  " 
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The  inclusion  of  Kq.  (28)  into  the  equation  set  for  a  family  of  samples  insures 
the  consistency  of  the  distribution  equations,  particularly  the  totals  equations, 
with  the  M  vs  Z  relation  that  is  presumed  for  the  family.  It  also  obviates  any 
need  to  determine  A  bv  least-squares  methods,  such  as  discussed  by  Smith  and 

i  11 

i  aro. 

When  we  first  began  using  the  truncated  equations  to  describe  the  approximate 
distribution  properties  of  hvdrometeors  along  trajectory  paths  of  re-entry  vehi¬ 
cles,  wc  had  virtually  no  information  about  the  nature  of  the  possible  variability 
of  the  AD  quantity  that  appears  in  the  equations  for  the  truncation  ratios.  Thus, 
we  made  the  obvious  first  assumption  that  this  quantity  might  have  a  constant 
value,  that  is, 

AD|n  -  C  .  (29) 

From  surface  disdrometer  data  for  rain  that  were  acquired  at  Wallops  Island, 
Virginia  (Plank),  “  we  deduced  that  the  value  of  the  constant  was  about  7.  5.  We 
also  recognized  that,  for  rain  containing  drops  of  the  breakup  size,  the  upper 
truncation  situation  would  no  longer  be  governed  by  Kq.  (29)  but  would  be  specified 
by 


D  -  breakup  diameter  s  5  turn  .  (30) 

In  the  fall  of  1981,  we  undertook  an  investigation  of  the  details  of  ADm  vari¬ 
ability  and  found  that  for  rain,  the  quantity  ADm  was  primarily  dependent  on  sam¬ 
pling  volume  and  was  secondarily  dependent  on  variations  of  liquid-water -content 
within  the  samples.  The  example  of  the  relatively  homogeneous  situation  of  rain 
that  was  aircraft  sampled  on  23  February  1977  (Figure  7)  depicts  the  way  that 
AD|(|  (solid  curve)  increased  with  sampling  interval  or  with  sampling  volume.  The 
shaded  envelope  on  either  side  of  the  curve  indicates  the  first  standard  deviation 
of  the  individual  ADm  values.  This  scatter  is  mostly  due  to  the  variations  of 
liquid-water-content  within  the  samples. 

Figure  7  indicates  that  the  ADm  values  increase  with  sampling  volume  and 
that  they  appear  to  attain  a  asymptotic  value  of  about  8.2  for  a  sampling  volume 

11.  Smith,  Jr.,  P.  L. ,  and  Laco,  C.  P.  (1978)  Techniques  for  fitting  size  distri¬ 

bution  functions  to  observed  particle  size  data.  Preprints  of  papers  pre¬ 
sented  at  the  18th  Conference  on  Radar  Meteorology.  Atlanta,  Georgia, 
28-31  March  1978,  pp.  129-133, 

12.  Plank,  V.G.  (1977)  Hydrometeor  data  and  analytical-theoretical  investiga¬ 

tions  pertaining  to  the  SAMS  Missile  Flights  of  the  1972-73  season  at 
Wallops  Island,  Virginia,  Environmental  Research  Papers  No,  603, 
AFGL/SAMS  Report  No.  5,  AFGL-TR -77 -01 49,  ADA051192,  239  pp. 
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somewhat  in  excess  of  10m  .  This  volume,  it  should  be  noted,  is  the  volume  of 


a  representative  atmospheric  sample  for  the  Alabama  rain  situation.  Any  sample 
of  appreciably  smaller  size  is  a  nonrepresentative  sample. 

Figure  8  compares  the  ADm  values  that  were  obtained  on  23  February  1977, 
with  those  from  two  other  rain  situations  sampled  by  the  MC-130R,  one  near  the 
Kwajalein  atoll  on  4  July  1978,  and  the  other  near  Hanscom  Air  Force  Base  on  1  5 
August  1979.  In  this  figure,  the  AD  ^  parameters  are  plotted  against  sampling 
volume.  The  latter  two  cases  were  of  much  shorter  duration  thus,  less  sampling 
volume. 
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Figure  8.  Fdot  of  AD  vs  Sampling  Volume  for  Three  Haiti 
Situations 

Although  these  latter  cases  display  considerably  more  variability  than  23 
February,  all  three  curves  show  small  AD^  values  at  small  sampling  volumes 
and  all  increase  markedly  up  to  a  volume  of  ~  6  m'^.  It  appears  that  an  asvmp- 

3 

totic  value  is  reached  at  —  10  m'  . 
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The  fin-ling  that  a  representative  atmospheric  sample  for  rain  is  of  the  order 
>f  1U  nr'  corresponds  with  the  comments  and  equations  of  Joss  and  Waldvogel^3 
in  that,  for  surface  rain  with  rates  from  1  to  10  mm  hr  \  about  ti  to  19  m3  sam¬ 
pling  volume  is  required  to  determine  radar  7.  values  from  size  distribution  data 

w  ith  a  0.9a  probability  of  90  pet  cent  accuracy.  Moreover,  Plank,  llerthel,  and 
1 1  3 

Delgado  “  have  shown  that  flight  durations  of  some  50  to  100  sec  (f>  to  10  m  vol¬ 
ume)  are  required  to  obtain  stable,  representative  size  distribution  information 
for  in. 

5.  COMPARISONS  OK  TRUNCATED  AND  NONTRUNCATF.D  DISTRIBUTIONS 
OK  EXPONENTIAL  TYPE 

To  illustrate  the  differences  between  truncated  exponential  distributions  and 
nontruncated  distributions,  we  hav  -  assumed  the  M  vs  /.  relation 

M  0.00314  7°-  :,7fi  g  m*3  ,  (ill 

which  corresponds  to  the ’’widespread  rain"  situations  of  Joss,  Thnms,  and 
Waldvogel13  and  which  wc  have  found  to  provide  good  descriptions  for  such  rains 
at  Wallops  Island  and  K/  ajalcin. 

From  this  assumed  M  vs  7.  relation  containing  a  total  liqui  i-w  atcr -content 
value  of  0.  2  g  m  ,  the  equation  set  described  pr  eviously  herein  was  solved  fot 
the  casi>  of  nontrunention.  In  such  case,  the  truncation  ratios  t  ^.  and  r^  and  \  y 
all  have  the  value  unity.  The  distribution  curves  for  number-  concent  ration  N,  for 
liquid-w-ater-eontent  M  and  for  radar  reflectivity  factor  Z  are  shown  in  Figure  9. 
They  are  the  solid  curves. 

For  comparison,  we  considered  two  other  situations.  First,  was  a  situation 
of  representative  atmospheric  sampling  for  which,  (in  the  Alabama  case)  AtJ  ‘ 
8.2,  r  =  0.090,  r^j  -  0.951'.,  and  r y  =  0.088.  The  equation  set  was  solved  for 
these  values  for  d  0.  2  mm  and  for  a  total  liquid-water-content  value  of  0.2 
g  m  ,  as  before.  The  distribution  curves  for  this  situation  are  also  shown  in 
Figure  9.  They  are  the  dashed  curves. 

We  next  considered  a  situation  of  typical  aircraft  sampling  with  PIUS  1-D 
instruments.  We  chose  a  5-sec  sampling  interval  that  corresponds  to  about  500  m 

13.  Joss,  J.,  and  Waldvogel,  A.  (1909)  Raindrop  sire  distribution  and  sampling 

size  distribution  and  sampling  size  errors,  J.  Atmos.  Sci.  20:500-509. 

14.  Joss,  J..  Thnms,  J.C.,  and  Waldvogd,  A.  (1908)  The  variation  of  raindrop 

size  distributions  at  Locarno,  Proc,  Intcrnatl.  Conf.  on  Cloud  Physics, 
Toronto,  Atru-r.  Meteorol.  Soc.  Boston,  p.  309. 
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Figure  9.  The  Distributed  Values  of  M,  and  Z  for  a  Nont runonted  Situation 
(Solid  Line)  and  Those  Truncated  with  DmA  -  4.5  (Dotted  Line)  and  Dm  A  3.  2 
(Dashed  I.ine).  In  all  cases.  M  0.2  (  m"3  and  7.  =  Dali 

of  flight  path  distance  at  sampling  airspeeds.  The  sampling  volume  for  this  inter¬ 
val  is  about  0.3  m  ,  which  is  considerably  smaller  than  a  representative  sam¬ 
pling.  With  reference  to  Figure  7,  the  D^  A  value  for  this  sampling  volume  is 

about  4.  5,  the  r., ,  r, .  ,  and  r_  are  about  0.  740,  0.  (>39,  and  0.  192.  rpspeet  i  velv. 

.\  A  l 

The  equation  set  was  solved  for  these  values  and  the  results  are  presented  in 
Figure  9  (dotted  curves). 

It  is  seen  that  the  differences  between  the  three  solutions  are  considerable. 

There  are  no  differences  in  the  M  values,  since  these  were  assumed  equal.  There 

are  also  no  differences  in  the  Z  values,  because  of  Lq.  (31).  However,  there  are 

maior  differences  in  the  A.  D  ,  and  values  that  are  listed  in  Table  1. 

m  T 
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Table5  1.  Comparison  of  the  A,  Dm,  and  Ny  from  a  Non- 
truncated  Distribution  and  Two  Truncated  Situations  with 
Different  Sampling  Volumes.  (M  -  0,2  f  m"3,  Z  =  1355.  6 
mm^  ;n  an  cases) 


*-l 

mm 

I'm 

mm 

N-r  3 

No.  m 

N’ont  runcat  »d 

AD  --  * 

3.23 

s: 

2150 

m 

Sampling  Volume  10  m'^ 

AD  =  8.  2 

2.  90 

2.  70 

1130 

m 

Sampling  Volume  0.3  m' 

AD  4.5 

|  1.  88 

2.  13 

550 

HI 

6.  SIMMARY 

In  this  report  we  have  demonstrated  the  conformity  of  rain  distributions  to 
exponential  shape  with  the  conformity  being  dependent  upon  sampling  volume.  We 
have  also  presented  a  set  of  equations  or  model  based  on  the  exponential  distribu¬ 
tion  function  that  may  be  utilized  for  estimating  the  parameters  of  rain  distribu¬ 
tions.  From  our  investigation  of  aircraft-acquired  measurements,  we  have  de¬ 
veloped  a  new  parameterization  i  ntity  (AD  )  that  appears  to  be  predictable  based 
on  sampling  volume. 

In  any  scenario  where  distributions  are  to  be  estimated,  some  assumption 
or  measurement  has  to  be  made  in  order  to  define  the  amount  of  precipitable 
water  in  a  given  volume  of  air.  We  chose  to  use  M  and  Z  as  the  definable  quan¬ 
tities  in  our  model  since  M  should  be  a  common  output  from  any  large-scale 
prediction  model  and  Z  can  be  defined  by  using  existing  M  vs  Z  relations  for  rain. 
Conversely,  definition  of  Z  can  be  made  by  remote  sensing  radar  and  evaluation 
of  1U  then  determined  through  these  same  relationships.  Introduction  of  M  and  Z 
into  this  model,  along  with  a  specified  sampling  volume  that  dictates  the  AD 
value  (Figure  8),  will  produce  reasonable  rain  number  and  size  information. 

The  number  and  sizes  of  raindrops  vary  considerably  both  within  a  specific 
situation  ami  between  situations.  As  examples,  the  variations  in  distribution 
parameters  for  the  three  cases  considered  in  this  study  are  presented  in  Appen¬ 
dix  A. 

This  model,  or  set  of  equations,  used  in  conjunction  with  actual  or  model 
output  values  of  M  and/or  Z  provides  raindrop  size  distributions.  These  distribu¬ 
tions  are  needed  for  theoretical  calculations  of  electromagnetic  rain  attenuation 
and  for  studies  of  rain  erosion  on  hypersonic  vehicles. 
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Appendix  A 

Variabilities  in  Distribution  Parameters  for  Three  Situations 

The  tables  (Tables  A1  through  A18)  presented  in  this  section  demonstrate  the 
variabilities  in  Dm,  ftl,  Z,  A,  and  AD|n  that  existed  in  the  three  situations 

sampled  in  »his  investigation.  They  are  included  to  provide  an  insight  as  to  the 
possible  scatter  that  may  exist  in  these  parameters  under  differing  conditions. 

XXp  and  Dm  values  are  directly  from  the  PMS  1-D  measurements.  M  and  Z 
are  calculated  from  these  measurements  as  explained  in  Section  3.  The  A  values 
were  derived  using  Kq.  (20)  and,  of  course,  the  A  Dpl  quantities  are  the  products 
of  these  two  parameters. 

There  is  a  table  for  each  of  the  parameters  listed  for  each  of  the  three  situa¬ 
tions.  The  tables  arc  divided  into  subsets  according  to  averaging  time  intervals. 
The  left  side  of  the  heading  on  each  subset  lists  the  averaging  time  in  seconds,  the 
number  of  averaged  points  in  that  particular  analysis,  the  number  of  points  where 
there  were  no  data  on  the  1-D  record  or  where  Kq.  (20)  could  not  be  solved  using 
the  particular  M,  Z,  and  1)  values,  the  total  number  of  1-sec  1-D  points  sam¬ 
pled,  ami  the  approximate  sampling  volume.  The  right  side  lists  the  mean,  min¬ 
imum,  ami  maximum  values  of  that  particular  number  oT  averaged  points  and  the 
standard  deviation  of  the  set. 

The  tabulation  below  each  heading  ranks  the  number  and  percentage  of  aver¬ 
aged  points  in  percentile  classes  relative  to  the  mean  value.  For  example,  in  the 
first  subset  of  Table  A2  the  averaging  time  was  a  see,  there  were  2Ht!  averaged 


points  and  since  all  the  1-sec  records  contained  data,  the  total  number  of  points 
sampled  were  1440.  Of  the  averaged  points,  44  or  ~  15  percent  of  the  288  were 
within  10  percent  of  the  mean  value.  Four  points  or  ~  1  percent  of  the  averaged 
points  were  larger  than  100  percent  from  the  mean. 
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Table  A 1.  Variability  in  on  23  February  1977 
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Table  A3.  Variability  in  M  on  23  February  1977 
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Variability  in  A  on  23  February  1977 
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Table  Afi.  Variability  in  AD  on  23  February  1977 
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Table  A  9.  Variability  in  M  on  4  July  1978 
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Table  A13.  Variability  in  N-,  on  15  August  1979 
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(Nnj,  =  No.  m  ,  sampling  volume  =  m  ) 
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Table  A14.  Variability  in  Dm  on  15  August  1979 
(Dm  =  mm,  sampling  volume  =  m^) 
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Table  A 15.  Variability  in  M  on  15  August  1979 
(M  -  sampling  volume  =  m‘*) 
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